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The hadroproduction of lepton pairs with mass Q and finite transverse momentum Qt is described in perturbative QCD by 
the same partonic subprocesses as prompt photon production. We demonstrate that, like prompt photon production, lepton 
pair production is dominated by quark-gluon scattering in the region Qt > Q/2. This feature leads to sensitivity to the 
gluon density in kinematical regimes accessible in collider and fixed target experiments, and it provides a new independent 
method for constraining the gluon density. 



1. INTRODUCTION 

The production of lepton pairs in hadron collisions 
/iift.2 7*^5 7* ~* II proceeds through an intermedi- 
ate virtual photon via qq — > 7*, and the subsequent 
leptonic decay of the virtual photon. Traditionally, in- 
terest in this Drell-Yan process has concentrated on 
lepton pairs with large mass Q which justifies the ap- 
plication of perturbative QCD and allows for the ex- 
traction of the antiquark density in hadrons [ ^ . 

Prompt photon production hih2 ■^X can be cal- 
culated in perturbative QCD if the transverse mo- 
mentum Qt of the photon is sufficiently large. Be- 
cause the quark-gluon Compton subprocess is domi- 
nant, gq — !• jX, this reaction provides essential in- 
formation on the gluon density in the proton at large 
X [ ^. Unfortunately, the analysis suffers from frag- 
mentation, isolation, and intrinsic transverse momen- 
tum uncertainties. Alternatively, the gluon density can 
be constrained from the production of jets with large 
transverse momentum at hadron colliders [ ^ , but the 
information from different experiments and colliders is 
ambiguous. 

In this paper we demonstrate that, like prompt pho- 
ton production, lepton pair production is dominated by 
quark-gluon scattering in the region Qt > Q/2. This 
realization means that new independent constraints on 
the gluon density may be derived from Drell-Yan data 
in kinematical regimes that are accessible in collider 
and fixed target experiments but without the theo- 
retical and experimental uncertainties present in the 
prompt photon case. 

In Sec. B, we review the relationship between vir- 



*Supported by the U.S. Department of Energy, Division of High 
Energy Physics, under Contract W-31-109-ENG-38. 
t Supported by Bundesministerium fiir Bildung und Forschung 
under Contract 05 HT9GUA 3, by Deutsche Forschungsgemein- 
schaft under Contract KL 1266/1-1, and by the European Com- 
mission under Contract ERBFMRXCT980194. 



tual and real photon production in hadron collisions in 
next-to-leading order QCD. In Sec. ^ we present our 
numerical results, and Sec. ^ is a summary. 

2. NEXT-TO-LEADING ORDER QCD FOR- 
MALISM 

In leading order (LO) QCD, two partonic subpro- 
cesses contribute to the production of virtual and real 
photons with non-zero transverse momentum: qq — > 



7(*)g and qg 7^*^?. The cross section for lepton 
pair production is related to the cross section for vir- 
tual photon production through the leptonic branch- 
ing ratio of the virtual photon a/ {^ttQ"^) . The virtual 
photon cross section reduces to the real photon cross 
section in the limit Q^ 0. 

The next-to-leading order (NLO) QCD corrections 
arise from virtual one-loop diagrams interfering with 
the LO diagrams and from real emission diagrams. 
At this order 2-^3 partonic processes with incident 
gluon pairs (gg) , quark pairs (qq) , and non-factorizable 
quark-antiquark (552) processes contribute also. Sin- 
gular contributions are regulated in n==4-2e dimensions 
and removed through MS renormalization, factoriza- 
tion, or cancellation between virtual and real contri- 
butions. An important difference between virtual and 
real photon production arises when a quark emits a 
collinear photon. Whereas the coUinear emission of a 
real photon leads to a 1/e singularity that has to be 
factored into a fragmentation function, the collinear 
emission of a virtual photon yields a finite logarithmic 
contribution since it is regulated naturally by the pho- 
ton virtuality Q. In the limit Q^ the NLO virtual 
photon cross section reduces to the real photon cross 
section if this logarithm is replaced by a l/e pole. A 
more detailed discussion can be found in [ Q . 

The situation is completely analogous to hard photo- 
production where the photon participates in the scat- 
tering in the initial state instead of the final state. For 
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real photons, one encounters an initial-state singular- 
ity that is factored into a photon structure function. 
For virtual photons, this singularity is replaced by a 
logarithmic dependence on the photon virtuality Q [ 

A remark is in order concerning the interval in Qt 
in which our analysis is appropriate. In general, in 
two-scale situations, a series of logarithmic contribu- 
tions will arise with terms of the type o^Io^XQ/Qt)- 
Thus, if either Qt » Q or Qt « Q, resumma- 
tions of this series must be considered. For practical 
reasons, such as event rate, we do not venture into 
the domain Qt » Q, and our fixed-order calcula- 
tion should be adequate. On the other hand, the cross 
section is large in the region Qt « Q. In previous pa- 
pers [ Q , we compared our cross sections with available 
fixed-target and collider data on massive lepton-pair 
production, and we were able to establish that fixed- 
order perturbative calculations, without resummation, 
should be reliable for Qt > Q/2. At smaller values of 
QT^ non-perturbative and matching complications in- 
troduce some level of phenomenological ambiguity. For 
the goal we have in mind, viz., contraints on the gluon 
density, it would appear best to restrict attention to 
the region Qt > Q/2, but below Qt » Q. 

3. PREDICTED CROSS SECTIONS 

In this section we present numerical results for the 
production of lepton pairs in pp collisions at the Teva- 
tron with center-of mass energy VS* — 1.8 and 2.0 
TeV. We analyze the invariant cross section Ecfia/dp^ 
averaged over the rapidity interval -1.0 < y < 1.0. 
We integrate the cross section over various intervals 
of pair-mass Q and plot it as a function of the trans- 
verse momentum Qt- Our predictions are based on 
a NLO QCD calculation [ ||] and are evaluated in the 
MS renormalization scheme. The renormalization and 
factorization scales are set to /i = /i/ = \/Q^ + Qt- 
If not stated otherwise, we use the CTEQ4M parton 
distributions [ Q and the corresponding value of A in 
the two-loop expression of as with four fiavors (five if 
/i > rrib). The Drell-Yan factor a/ (SttQ^) for the decay 
of the virtual photon into a lepton pair is included in 
all numerical results. 

In Fig. |l| we display the NLO QCD cross section for 
lepton pair production at the Tevatron at =1.8 
TeV as a function of Qt for four regions of Q. The 
regions of Q have been chosen to avoid resonances, 
i.e. between 2 GeV and the J/ip resonance, between 
the J/4> and the T resonances, above the T's, and a 
high mass region. The cross section falls both with 
the mass of the lepton pair Q and, more steeply, with 
its transverse momentum Qt- No data are available 




Figure 1. Invariant cross section Ecfia/dp^ as a func- 
tion of Qt for pp j*X at ^/S = 1.8 TeV in non- 
resonance regions of Q. The cross section falls with 
the mass of the lepton pair Q and, more steeply, with 
its transverse momentum Qt- 



yet from the CDF and DO experiments. However, 
prompt photon production data exist to Qt — 100 
GeV, where the cross section is about 10~^ pb/GeV^. 
It should be possible to analyze Run I data for lep- 
ton pair production to at least Qt — 30 GeV where 
one can probe the parton densities in the proton up to 
XT = 2Qt/VS ~ 0.03. The UAl collaboration mea- 
sured the transverse momentum distribution of lepton 
pairs at ^/S — 630 GeV up to xt — 0.13 [||], and their 
data agree well with our theoretical results [ 

The fractional contributions from the qg and qq sub- 
processes through NLO are shown in Fig. ||. It is evi- 
dent that the qg subprocess is the most important sub- 
process as long as Qt > Q/2. The dominance of the 
qg subprocess diminishes somewhat with Q, dropping 
from over 80 % for the lowest values of Q to about 70 % 
at its maximum for Q ~ 30 GeV. In addition, for very 
large Qt, the significant luminosity associated with the 
valence dominated q density in pp reactions begins to 
raise the fraction of the cross section attributed to the 
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Figure 2. Contributions from the partonic sub- 
processes qg and qq to the invariant cross section 
E(Pa/dp^ as a function of Qt for pp at \/S 

~ 1.8 TeV. The qg channel dominates in the region 
Qt > Q/2. 



Figure 3. Invariant cross section Ecfia/dp^ as a func- 
tion of Qt for pp ^ and two different center-of- 
mass energies of the Tevatron (Run 1: = 1.8 TeV, 
Run 2: V5 = 2.0 TeV). The cross section for Run 2 is 
5 to 20 % larger, depending on Qt- 



qq subprocesses. Subprocesses other than those initi- 
ated by the qq and qg initial channels are of negligible 
import. 

We update the Tevatron center-of-mass energy to 
Run II conditions {\fS = 2.0 TeV) and use the latest 
global fit by the CTEQ coUaboration (5M). Figure | 
demonstrates that the larger center-of-mass energy in- 
creases the invariant cross section for the production of 
lepton pairs with mass 5 GeV < Q < 6 GeV by 5 % at 
low Qt^I GeV and 20 % at high Qt ^ 100 GeV. In 
addition, the expected luminosity for Run II of 2 fb~^ 
should make the cross section accessible to Qt — 100 
GeV or xt — 0.1. This extension would constrain the 
gluon density in the same regions as prompt photon 
production in Run I. 

Next we present a study of the sensitivity of collider 
and fixed target experiments to the gluon density in 
the proton. The full uncertainty in the gluon density 
is not known. Here we estimate this uncertainty from 
the variation of different recent parametrizations. We 



choose the latest global fit by the CTEQ collaboration 
(5M) as our point of reference [ || and compare results 
to those based on their preceding analysis (4M[0]) and 
on a fit with a higher gluon density (5IIJ) intended 
to describe the CDF and DO jet data at large trans- 
verse momentum. We also compare to results based on 
global fits by MRST [ |^ , who provide three different 
sets with a central, higher, and lower gluon density, 
and to GRV98 [ ||]. 

In Fig. U we plot the cross section for lepton pairs 
with mass between the J /ij) and T resonances at Run II 
of the Tevatron in the region between Qt = 10 and 30 
GeV {xT = 0.01 . . . 0.03). For the CTEQ parametriza- 
tions we find that the cross section increases from 4M 
to 5M by 2.5 % {Qt = 30 GeV) to 5 % {Qt = 10 

^In this set a purely perturbative generation of heavy flavors 
(charm and bottom) is assumed. Since we are working in a 
massless approach, we resort to the GRV92 parametrization for 
the charm contribution [|l^ and assume the bottom contribution 
to be neghgible. 
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Figure 4. Invariant cross section ESa/dp^ as a func- 
tion of Qt for pp -f*X at^/S^ 2.0 TeV in the re- 
gion between the J/ip and T resonances. The largest 
differences from CTEQ5M are obtained with GRV98 
at low Qt (minus 10 %) and with MRST(gt) at large 
Qt (minus 7 %). 



GeV) and from 5M to 5HJ by 1 % in the whole Qr- 
range. The largest differences from CTEQ5M are ob- 
tained with GRV98 at low Qt (minus 10 %) and with 
MRST(gT) at large Qt (minus 7%). 

The theoretical uncertainty in the cross section can 
be estimated by varying the renormalization and fac- 
torization scale fi = fif around the central value 
\/ Q'^ + Qt- Figure || shows this variation for pp 
j*X at VS = 2.0 TeV in the region between the 
J/ip and T resonances. In the interval 0.5 < 
m/ y/Q"^ + Qt < dependence of the cross section 

on the scale — fif drops from ±15% (LO) to the small 
value ±2.5% (NLO). The i^T-factor ratio (NLO/LO) is 
approximately 2, as one might expect naively. 

A similar analysis for Fermilab's fixed target experi- 
ment E772 [ is shown in Fig. ^. In this experiment, 
a deuterium target is bombarded with a proton beam 
of momentum piab = 800 GeV, i.e. yfS = 38.8 GeV. 
The cross section is averaged over the scaled longitu- 
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Figure 5. Invariant cross section Ecfia/dp^ as a func- 
tion of the renormalization and factorization scale fi = 
jjLf for pp ^ 7*X at =2.0 TeV in the region be- 
tween the J/ tp and T resonances and Qt = 5.5 GeV. In 
the interval 0.5 < fJ,/ y/Q'^ + Qt < 2 the dependence of 
the cross section on the scale ^ — drops from ±15% 
(LO) to ±2.5% (NLO). The X-Factor (NLO/LO) is 
approximately 2. 



dinal momentum interval 0.1 < xp < 0.3. In fixed 
target experiments one probes substantially larger re- 
gions of XT than in collider experiments. Therefore one 
expects greater sensitivity to the gluon distribution in 
the proton. We find that use of CTEQ5HJ increases 
the cross section by 7 % (26 %) w.r.t. CTEQ5M at 
Qt = 3 GeV {Qt = 6 GeV) and by 134 % at Qt = 10 
GeV. With MRST(g|) the cross section drops relative 
to the CTEQ5M-based values by 17 %, 40 %, and 59 
% for these three choices of Qt- 

Figure shows the variation of the fixed target cross 
section on the renormalization and factorization scale 
fi — ^f. In the interval 0.5 < m/VQ^ + Qt 
the dependence decreases from ±49% (LO) to ±37% 
(NLO). An optimal scale choice might be /i = /i/ = 
\/Q^ + Qt/^1 where the points of Minimal Sensitiv- 
ity (maximum of NLO) and of Fastest Apparent Con- 
vergence (LO=NLO) nearly coincide. At /x = /i/ = 
+ Qt' if -factor ratio is 2.6. The NLO cross 
section turns negative at the lowest scale shown /i = 
f^f — + Qt/8 — 1 GeV, a value too low to guar- 
antee perturbative stability. 
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Figure 6. Invariant cross section EcPa/dp'^ as a func- 
tion of Qt for pN 7*X at pu^ = 800 GeV. The 
cross section is highly sensitive to the gluon distribu- 
tion in the proton in regions of xt where it is poorly 
constrained in current analyses. 



4. SUMMARY 

The production of Drell-Yan pairs with low mass 
and large transverse momentum is dominated by gluon 
initiated subprocesses. In contrast to prompt pho- 
ton production, uncertainties from fragmentation, iso- 
lation, and intrinsic transverse momentum are absent. 
The hadroproduction of low mass lepton pairs is there- 
fore an advantageous source of information on the 
parametrization and size of the gluon density. The in- 
crease in luminosity of Run II increases the accessible 
region of xt from 0.03 to 0.1. The theoretical uncer- 
tainty has been estimated from the scale dependence 
of the cross sections and found to be very small for 
collider experiments. 

Acknowledgment 

It is a pleasure to thank L. E. Gordon for his collab- 
oration. 



X 10 

0.3 p- 

0.25 - 

0.2 - 



pN -> y X at pj^i, = 800 GeV 



5 GeV < Q < 6 GeV: 
Q.^ = 5.5 GeV _! 
0.1<Xp<0.3 : 

NLO : 

LO - 




ivV(q'+Qt ) 



Figure 7. Invariant cross section Ed^a/dp^ as a func- 
tion of the renormalization and factorization scale /i = 
/i/ for pN 7*X at puh = 800 GeV. In the interval 
0.5 < [ij \/Q'^ + Q\ < 2 the dependence of the cross 
section on the scale ^ drops from ±49% (LO) to ±37% 
(NLO). 
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